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We use the continuum shell model approach to explore the resonance width distribution in unstable 
many-body systems. The single-particle nature of a decay, the few-body character of the interaction 
Hamiltonian, and collectivity that emerges in non-stationary systems due to the coupling to the 
continuum of reaction states are discussed. Correlations between structures of the parent and 
daughter nuclear systems in the common Fock space are found to result in deviations of decay width 
statistics from the Porter-Thomas distribution. 
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The Porter-Thomas distribution (PTD) [1] of transi- 
tion strengths is a central aspect of complex systems. 
This statistical law was noted by many authors [2] to be 
valid more generally than other predictions of the Ran- 
dom Matrix Theory from which it originates. The PTD 
emerges under the assumption that the relative orienta- 
tion of the two states involved in the overlap describing a 
transition covers the ri-dimensional sphere in the Hilbert 
space uniformly. While this is true by definition for the 
Gaussian Orthogonal Ensemble (GOE), the validity of 
the PTD extends much farther as it constitutes the cen- 
tral limit theorem (CLT). Being a sum of a large number 
of uncorrelated components, the transitional amplitude is 
indeed expected to have a Gaussian distribution. There 
is a large volume of work on this subject; see reviews \2r 
Q and references therein. Generally, there is a consensus 
among authors that while the specifics of an ensemble 
and the physics of transitions do matter for certain ob- 
servables, any deviations from the PTD are quickly de- 
feated by even small stochastic components due to the 
robust nature of the CLT; see for example Ref. @- Any 
claims to the contrary, either experimental @ or theoret- 
ical 0, [Tol | , have always ignited debates and discussions 

Here we will not consider the data handling procedures, 
which on many occasions, have been deemed to be the 
most likely reasons for the deviations observed experi- 
mentally [6j. Instead, we focus on the possible reasons 
for high correlation between the transitioning states from 
the theoretical perspective. We analyze the feasible sce- 
narios with the help of the continuum shell model ap- 
proach [no which is one of the most equipped meth- 
ods to address the structure-reaction physics of interest 
microscopically. In a unified picture we review the super- 
radiance (SR) effects [fj, ll4( , the wave-function localiza- 
tion effects in the two-, three-, and four-body ensembles 
(lTjj . the role of rotational symmetry, and other parent- 
daughter structural correlations that emerge in a decay 

The dynamics of an unstable many-body system pro- 
jected onto the intrinsic space spanned by the bound 
(shell-model) states is generated by the effective, energy- 



dependent Hamiltonian [H Oil [13 

H = H-i ]T p c \c)(c 

c(opcn) 



(1) 



Here H is the Hermitian part that is identified with the 
traditional shell model Hamiltonian. The second, imag- 
inary term reflects the irreversible decays into the con- 
tinuum of states excluded by the Feshbach projection. 
This factorized operator contains the kinematic penetra- 
bility factor p c and the set of channel vectors |c). For 
simplicity we omit the angular momentum, isospin, and 
other labelings; detailed notations are found in Ref. [T^ |. 
The problem is non-stationary; the Hamiltonian (JXJ) is 
understood as a component of the propagator and is de- 
pendent on the scattering energy. The Hermitian com- 
ponent includes the coupling to the continuum of reac- 
tion states via virtual excitations; the penetrability also 
depends on energy through the kinematics of the decay 
process. Away from thresholds the energy dependence is 
smooth, and its exact form mainly pertains questions of 
the experimental data analysis. We ignore this depen- 
dence here, and further assume that the penetrability is 
the same for all channels. The eigenvalues of the effec- 
tive Hamiltonian (Q]) are complex, £ = E — iT/2, and 
represent the poles of the scattering matrix in the com- 
plex energy plane. These complex energies are associated 
with resonances and their widths. 

Let us first consider weak decays, for which the imag- 
inary component in ([T]) can be treated perturbatively. 
In this case the shell model eigenstate \I) defined by 
H\I) = Ej\I) is not modified by the decay instability, 
and the corresponding decay width is 

r/ = 2p7/, where 7j = ^ |(/|c)| 2 

c(opcn) 

is the reduced width. The PTD of reduced widths 
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emerges under the uniform Hilbert space coverage as- 
sumption for | J). Here v is the dimension spanned by the 
channel vectors, and 7 is the average reduced width. For 
the orthogonal and normalized channels 7 = v /O. 



2 



In this work we examine situations with only one open 
channel. The strength of the continuum coupling is de- 
fined via the average decay width relative to the level 
spacing. Here we express this coupling using the param- 
eter k — p/X, where A 2 = f2 _1 Tr (-ff 2 ) is the variance of 
the density of states distribution of H. 

In Fig. [JJwe consider an example of GOE+SR, where 
H in (dJ is represented by the GOE. Here the PTD is 
reproduced numerically in the limit p — > 0. The den- 
sity of states in this limit has a semicircular distribution 
bound by a radius 2A. The imaginary component in ([1]) 
is factorized, which is pertinent to the unitarity of the 
scattering matrix. The non-Hermitian component, when 
large, gives rise to the collectivity often referred to as 
superradiance (SR). A similar collectivity due to the fac- 
torized Hermitian interaction describes giant resonances. 
The resultin g d eformed random ensembles are discussed 
in Refs. (2I Tl8L Ojj|. As coupling to the continuum in- 
creases, and the average decay width relative to the level 
spacing k becomes large, the resonances start to overlap, 
thus reorienting the intrinsic structure. This could be 
hypothesized to result in the PTD being violated Q . An 
in-depth examination, however, shows that the SR mech- 
anism alone is unlikely to cause a significant change to the 
PTD. Indeed, for k <C 1, the PTD simply follows from 
the definition of GOE. For k ^> 1, in full mathematical 
equivalence to the deformed ensembles, the Hilbert space 
is separated into the SR channel space, which is one- 
dimensional here with a single eigenstate £ — —ip, and 
the orthogonal statistical (compound resonance) space of 
dimension f2 — 1. Because H is orthogonally invariant the 
reduced-space dynamics is represented by the GOE. With 
the perturbation theory built in this limit one finds that 
the reduced widths for the compound resonances follow 
the PTD with 7 = l/(f}/t 2 ). For large Q the single SR 
state with a reduced width jsr — 1 — K ~ 2 nas 110 effect on 
the PTD. As seen in Fig. [TJ numerical study confirms the 
PTD for both small and large values of k. Moreover, the 
slight deviation from the PTD for couplings k between 
around 0.4 and 2 is due to a small fraction of very broad 
states with 7 > IO7. This localized effect is shown in the 
inset of Fig. [TJ This subset of an exceptionally broad 
states is difficult to identify experimentally. Investiga- 
tions [2j, [lg, ll9j of deformed ensembles further confirmed 
a good agreement of the decay width distribution with 
the PTD and of the density of states distribution with 
the semicircular shape. 

Experience with the realistic nuclear structure and 
some theoretical arguments suggest that the effec- 
tive Hamiltonian involves only few-nucleon interactions, 
thus, the two-body random ensembles (TBRE) appear to 
be more appropriate. Many features of these ensembles 
are different from those of GOE, nevertheless the eigen- 
vectors form a uniform coverage of the Hilbert space. 
Numerical studies confirm that, in agreement with the 
CLT, this leads to the PTD of transition strengths to- 
ward an uncorrelated channel vector [7]. This logic, how- 
ever, does not take into account the correlations that ex- 
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Figure 1: (Color online) The normalized distribution of prob- 
abilities of the reduced widths for the GOE+SR ensemble 
with n = 10 4 . The curves for different continuum couplings 
k = 0.5, 1, 2, and 5 are compared with the PTD. The cor- 
responding average widths Qrf = 0.76, 0.40, 0.20, and 0.05, 
are computed after the very broad state(s) are removed from 
consideration. For « > 2 this amounts to exclusion of a single 
superradiant state in which case = k™ 2 . With the ex- 
ception of very broad states, there is perfect agreement with 
PTD, the differences being magnified in the inset where the 
region of 7 > 47 is shown. The largest deviation is observed 
for the transitional coupling strength k = 1. For k = 5 the 
distribution is already indistinguishable from the PTD. 

ist in the variable particle-number Fock space. Strong 
parent-daughter correlations emerge due to the micro- 
scopic physics of decay. Indeed, in the single-particle 
reaction processes all nucleons, except for one, are spec- 
tators, and the decay channels for the N— particle system 
c; N) are built from the (N — l)-particle eigenstates of 
the daughter nucleus \F; N — l) that follow from the same 
two-body Hamiltonian. Thus, \c;N) = {a<\F;N - 1)} , 
where a' is a single-particle creation operator correspond- 
ing to the decaying nucleon, and brackets {...} indicate 
normalization to unity and appropriate symmetry cou- 
pling. The correlation between eigenstates and channels 
leads to the violation of PTD. Single-particle removal 
amplitudes are related to the independent-particle basis, 
where departure from the PTD has been demonstrated 
in the past |10( . 

In Fig. [5]the EGOE+SR ensemble is considered, where 
H in Eq. JTJ) is given by an embedded two-body GOE 
(EGOE), and by definition does not include any sym- 
metries. The previous example shows that it is im- 
portant to draw attention to the statistics of the nar- 
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Figure 2: (Color online) The normalized distribution of prob- 
abilities of decay amplitudes in the EGOE+SR ensemble. The 
decay of a 7-particle system to the ground state of a 6-particle 
system is considered. There are 16 single-particle states, thus 
Q — 11440. The curves for different continuum couplings 
K = 0, 1, and 2 are compared with the Gaussian and Bessel 
distributions. The curves are normalized so that the average 
widths do not include the single SR state, ffy = 1, 0.19, and 
0.06 for k = 0, 1 and 2 respectively. The region of very nar- 
row widths is shown in the inset using a log-log scale. While 
all observed distributions for very narrow states seem to ap- 
proach a constant, they are still not described by Gaussian 
distribution of different variances (or 7's). The inset includes 
two Gaussian curves with variances 1 and 0.02 as labeled. 



row states. In this work we do not present our fits of 
the observed distributions to the PTD @ treating v as 
a parameter. We found that attempts to do so could 
produce misleading results. The effective v is directly 
related to the normalization 7 which is disproportion- 
ately influenced by statistically unimportant collective 
state(s). In an attempt to avoid any potential misjudge- 
ments, here we show the distribution of absolute values 
of amplitudes x = y/j/j- Then the Gaussian distribu- 
tion Pg{x) = \j2jit cxp (— x 2 /2) that corresponds to the 
PTD is easy to separate from a distribution provided 
by the Bessel function Pb{x) — (2/it)Kq(x). In con- 
trast to the CLT, the latter distribution emerges when 
the transitional overlap is possible only due to a sin- 
gle component in the wave-function along some direc- 
tion in the Hilbert space given by a vector |1); so that 
(I\c) = (I|l}(l|c), where both <7|1) and (l|c) are dis- 
tributed normally (agree with PTD). 

From Fig. H we find that none of the EGOE+SR re- 
sults follow the PTD. In contrast to the Gaussian curve 



Pg(x) the distributions have sharp peaks at low ampli- 
tudes and an extended exponential tail. In the SR limit 
of large n the distributions appear to approach the one 
given by the Bessel function Pb(x). 

The rank of the force beyond the two-body interaction, 
and symmetries, such as rotational, may have additional 
influences. To examine this we consider a model where N 
identical fermions occupy a single- j level. This has been 
a pop ular model for exploring the properties of TBRE [3, 
Rl2dj| . The resemblance of the low- lying spectra to those 
observed in realistic nuclei is the most intriguing feature. 
For our demonstration we select j = 19/2 and discuss 
widths of the decay of many-body states in 9-particle 
systems. The final state is the ground state of the system 
with 8 nucleons. All states, in both parent and daughter 
nuclei, are eigenstates of the same Hamiltonian given by 
the n— body Random Ensemble (n-BRE) [2l|. In this 
work we restrict our consideration to the two-, three-, 
and four-body forces, n = 2, 3, and 4, respectively. We 
select only those realizations where the daughter system 
has ground state spin F — 0; and thus the channel spin 
is I = 19/2. The fractions of such realizations are 42%, 
64%, and 83% for n = 2,3, and 4 respectively. 

The resemblance between random ensembles with sym- 
metries and realistic nuclei extends to parentage rela- 
tions. The low-lying states in the odd-particle parent 
nucleus are predominantly of the single-particle nature. 
If F = for the even-particle core then the ground state 
of a system with an extra nucleon is likely to carry the 
single-particle quantum numbers j = I = 19/2. This 
is indeed observed, and the corresponding probabilities 
are 21%, 47%, and 37% for n — 2,3, and 4 respectively. 
The correlation between the parent and daughter ground 
states is demonstrated in Fig. [3] which shows the distri- 
bution of reduced widths for the decay from ground state 
to ground state when F = and I = 19/2. Both parent 
and daughter systems have correlated structures because 
they are eigenstates of the same Hamiltonian for different 
number of particles. In the distribution of spectroscopic 
factors this correlation is seen as a peak near 7 = 1. As 
the rank of interaction n becomes higher, more remote 
configurations can be admixed, which reduces the ground 
state to ground state transitional collectivity. For totally 
uncorrelated systems the transitional strength due to the 
CLT is expected to follow the PTD. 

Since most of the transitional strength is concentrated 
in a few states at the low end of the spectrum, here the 
PTD is not generally expected; however, there is also no 
agreement with this law for the widths of highly excited 
compound resonances. In Fig. 2] the distribution of re- 
duced decay amplitudes is shown for the two-, three-, and 
four-body random ensembles with rotational symmetry. 
The curves are quite close to each other and are similar to 
the EGOE result (« = in Fig. |2J|. All findings indicate 
violation of the PTD. 

This study is motivated by the long-standing debate in 
relation to the Porter-Thomas distribution and possibil- 
ity of its violation [7Hll| . The PTD is a robust prediction 
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Figure 3: (Color online) Distribution of the reduced decay 
widths (spectroscopic factors) for the ground state N — 9 
spin I = 19/2 parent decaying to the ground state N = 8 spin 
F — daughter system. The n-BRE of identical nucleons in 
j — 19/2 level is considered. Three curves correspond to two- 
, three-, and four-body (n = 2, 3, and 4) random ensembles 
with rotational symmetry. 
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justified by the central limit theorem; it is easily con- 
firmed for different random matrix ensembles. This, how- 
ever, is a purely structural approach that does not take 
into account the microscopic physics of reactions. To ad- 
dress this we use a continuum shell model approach where 
violations of the PTD may result from one or a combina- 
tion of the following: coherence in structure due to fac- 
torized nature of the effective Hamiltonian that reflects 
unitarity of the scattering matrix, the so called super- 
radiance mechanism; parent-daughter relation between 
decaying systems in the common Fock space; few-body 
low rank interaction forces; and significant variations in 
the energy dependence of the effective Hamiltonian. We 
examine all of these possibilities, with the exception of 
the last, which is to be discussed elsewhere. 

In agreement with the studies of deformed ensembles 
we find that the SR mechanism by itself does not lead 
to violations of the PTD. However, the distribution, un- 
ambiguously different from PTD, is observed in random 
ensembles with a particle decay and with few-body intrin- 
sic Hamiltonians. This remains true even when subsets 
of states are considered: restricted by a certain energy 
region, or by some reasonable limits on the value of the 
decay width itself. The latter allows for any collective or 
unmeasurably broad states to be excluded. The parent- 
daughter relation in the decay process appears to be cen- 
tral for this phenomenon, the results being only slightly 
influenced by SR, additional rotational symmetry or the 
rank of forces. 



Figure 4: (Color online) Same model as in Fig. [3] The dis- 
tribution of spectroscopic amplitudes for all 204 states with 
spin I — 19/2 in the parent 9-nucleon system is shown in the 
decay to the spin F — ground state of the daughter nucleus. 
The results are compared with the Gaussian and Bessel distri- 
butions. The curves do not change if states only in a certain 
energy region are considered. 
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